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ABSTRACT 

In  this  paper  we  derive  a  method  for  finding  small  amplitude  high 
frequency  solutions  to  hyperbolic  systems  of  quasilinear  partial  differential 
equations.  Our  solution  is  the  sum  of  two  parts:  (i)  a  superposition  of 
small  amplitude  high  frequency  waves;  (ii)  a  slowly  varying  'mean  solution'. 
Each  high  frequency  wave  displays  nonlinear  distortion  of  the  wave  profile  and 
shocks  may  form.  Shock  conditions  are  derived  for  conservative  systems. 
Different  high  frequency  waves  do  not  interact  provided  the  frequencies  and 
wave  numbers  of  two  waves  are  not  linearly  related  to  those  of  a  third.  The 
mean  solution  is  found  by  solving  a  linear  partial  differential  equation. 

This  method  generalizes  Whitham's  nonlinearization  technique  [9]  for  single 
waves,  to  problems  where  many  waves  are  present.  We  obtain  these  results  by 
generalizing  a  scheme  first  proposed  by  Choquet-Bruhat  [1]  which  employs  the 
method  of  multiple  scales . 
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SIGNIFICANCE  AND  EXPLANATION 


Problems  involving  high  frequency  waves  arise  in  many  applications.  For 
example  they  arise  in  optics  (light  waves  propagating  through  a  medium  of 
slowly  varying  refractive  index),  acoustics  (sound  waves  propagating  through 
the  atmosphere)  and  oceanography  (water  waves  propagating  over  a  gently 
sloping  ocean  bed).  If  the  waves  are  sufficiently  small  in  amplitude  they 
satisfy  linear  differential  equations  and  there  is  a  powerful  theory  for 
calculating  their  behaviour  called  geometrical  optics.  Larger  amplitude  waves 
may  satisfy  nonlinear  equations,  and  then  they  display  qualitatively  different 
effects  from  linear  waves.  The  most  significant  effect  of  nonlinearity  is 
that  waves  can  form  shocks,  surfaces  across  which  the  wave  amplitude  changes 
sharply.  The  most  familiar  example  of  a  shock  is  the  sonic  boom  generated  by 
an  aircraft  in  supersonic  flight.  In  this  paper  we  develop  a  generalization 
of  geometrical  optics  to  small  amplitude  nonlinear  waves  which  allows  us  to 
calculate  how  they  propagate,  and  to  predict  effects  such  as  the  formation  of 
shocks. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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WEAKLY  NONLINEAR  HIGH  FREQUENCY  WAVES 


John  Hunter 

1 .  Introduction 

In  this  paper  we  consider  waves  which  satisfy  nonlinear  hyperbolic 
partial  differential  equations.  Work  on  such  waves  has  been  centered  on  one¬ 
dimensional  problems.  The  problems  either  involve  only  two  independent 
variables,  usually  time  and  one  space  variable,  or  involve  only  a  single  wave 
and  thus  are  essentially  one-dimensional. 

Here  we  develop  a  method  for  problems  in  any  number  of  dimensions  and 
with  any  number  of  waves  present.  We  call  our  method  weakly  nonlinear 
geometrical  optics,  since  it  reduces  to  geometrical  optics  for  linear 
systems.  It  applies  to  small  amplitude,  high  frequency  waves  which  are 
solutions  of  quasi-linear  hyperbolic  partial  differential  equations. 

We  obtain  the  first  term  in  an  asymptotic  expansion  for  such  a  solution, 
which  is  the  sum  of  two  parts:  (i)  a  slowly  varying  mean  solution;  and  (ii) 
the  superposition  of  a  number  of  small  amplitude,  high  frequency  waves.  The 
profile  of  each  high  frequency  wave  in  the  superposition  is  distorted  by  the 
nonlinear  self-interaction  of  the  wave.  This  can  cause  the  solution  for  the 
wave  to  become  multivalued.  To  obtain  a  single-valued  solution  we  must 
introduce  discontinuities,  and  for  conservative  systems  we  derive  an  equal 
area  rule  that  allows  us  to  fit  shocks  into  the  wave  solution. 

We  show  that  there  is  no  interaction  between  different  waves  in  the 
superposition  provided  that  (a)  the  mean  solution  is  correctly  chosen;  and  (b) 
a  certain  resonance  condition  never  holds.  The  mean  solution  must  be  chosen 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 .  This 
material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  No.  MCS-7927062,  Mod.  1. 


so  that  an  average  of  the  wave  profile  is  zero.  Thus  the  mean  solution  is  the 
average  of  the  first  term  in  the  asymptotic  expansion.  It  is  found  by  solving 
a  linear  system  of  partial  differential  equations. 

Resonance  can  occur  when  the  frequencies  and  wave  numbers  of  two  waves 
are  linearly  related  to  those  of  a  third.  This  threefold  resonance  is  a 
consequence  of  the  quadratically  nonlinear  interactions  between  the  waves,  and 
a  more  general  form  of  solution  than  that  considered  here  must  be  used  if 
there  is  resonance. 

The  method  by  which  we  obtain  these  results  is  an  extension  of  one 
proposed  by  Choquet-Bruhat  [1],  and  employs  the  method  of  multiple  scales.  It 
gives  exactly  the  same  result  for  a  single  wave  as  the  nonlinearization 
technique  proposed  by  Landau  [2]  and  Whitham  [3],  and  later  derived  by  Varley 
and  Cumberbatch  ( [4] ,  [5] )  using  the  method  of  relatively  undistorted  waves. 

The  result  that  two  small  amplitude  waves  do  not  interact  to  first  order 
in  the  wave  amplitude,  has  been  found  previously  in  one  dimension  by  Mortell 
and  Varley  [6]  and  Seymour  and  Mortell  [7] .  The  present  method  generalizes 
this  result  for  high  frequency  waves  to  any  number  of  dimensions. 

An  outline  of  the  contents  of  this  paper  are  as  follows.  In  section  2  we 
summarize  the  results  obtained  by  the  method  of  weakly  nonlinear  geometrical 
optics.  The  formal  derivation  of  the  equations  stated  in  section  2  for  the 
mean  solution  and  the  high  frequency  waves,  is  given  in  section  3.  The  equal 
area  rule  for  fitting  shocks  into  the  high  frequency  waves  is  derived  in 
section  4. 

In  the  next  two  sections  we  compare  our  results  with  previous  results. 

In  section  5  we  show  that  Whitham* s  nonlinearization  technique  is  a 
consequence  of  the  present  method.  This  demonstrates  the  agreement  between 
weakly  nonlinear  geometrical  optics,  when  applied  to  a  single  wave,  and  the 
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method  of  relatively  undistorted  waves,  when  expanded  for  small  amplitudes. 
Both  methods  reduce  to  the  nonlinearization  technique  in  the  common  case  of  a 
single  small  amplitude  wave. 

In  section  6  we  consider  time- dependent  plane  and  spherically  symmetric 
waves  in  gas  dynamics.  The  equations  governing  such  waves  involve  only  two 
independent  variables  and  therefore  can  be  treated  using  the  method  of 
characteristics.  We  verify  the  results  of  weakly  nonlinear  geometrical 
acoustics  in  this  case  by  showing  that  they  also  follow  from  the  method  of 
characteristics.  In  particular  we  see  that  two  waves  do  not  interact  if  the 
mean  solution  is  chosen  as  described  above. 

In  section  7  we  discuss  conditions  under  which  resonance  can  occur, 
giving  necessary  and  sufficient  conditions  for  it  to  happen.  This  gives  us 
the  conditions  under  which  a  solution  consisting  of  the  superposition  of  non¬ 
interacting  waves  is  applicable. 

In  section  8  we  consider  systems  with  multiple  characteristics  and  show 
how  the  results  above  are  modified,  assuming  that  the  multiplicity  of 
characteristics  does  not  change. 

Finally  in  the  appendix  we  apply  the  method  to  the  gas  dynamics 
equations.  The  result  is  a  theory  of  weakly  nonlinear  geometrical 
acoustics.  Applications  of  the  theory  will  be  described  in  a  future  paper. 

I  wish  to  thank  Professor  J.  B.  Keller  for  suggesting  this  problem  to  me 
and  for  many  helpful  discussions  during  my  work  on  it.  This  work  was 
performed  while  the  author  was  at  Stanford  University  and  I  am  happy  to 
acknowledge  their  support. 


2.  summary  of  Weakly  Nonlinear  Geometrical  Optics 


We  consider  a  strictly  hyperbolic  system  of  quasi-linear  first  order 
partial  differential  equations 

(2.1)  £  A^^(x,u)u  +  B( x , u )  “0 

i=1  Xi 

In  (2.1)  x  =  (x1,...,xn)  is  an  n-vector,  B  and  u  =  (u.,,...,^)  are  N- 
vectors  and  each  A^^  is  an  N  *  N  matrix.  Usually  one  of  the  independent 
variables  x^  is  time.  We  assume  that  A^^  and  B  are  continuously 
differentiable  with  respect  to  u. 

The  weakly  nonlinear  geometrical  optics  solution  to  (2.1)  is 
m 

(2.2)  u  -  u(x,e)  +  e  £  a*-* *  (x,$^ * (x)/e)R(^ *  (x)  +  0(e2)  (as  e  o)  . 

j-1 

The  derivation  of  (2.2)  is  given  in  section  3.  Here  we  explain  what  the 
various  functions  appearing  in  (2.2)  are  and  how  they  are  determined. 

Equation  (2.2)  represents  u  as  the  sume  of  a  mean  solution  u  and  a 
superposition  of  m  small  amplitude  high  frequency  waves  The 

small  parameter  e  is  introduced  into  the  problem  through  initial  and 
boundary  conditions.  It  is  the  ratio  of  a  typical  dimensioned  wave  amplitude 
relative  to  a  parameter  with  the  same  dimensions  appearing  in  the  problem,  it 
is  also  the  ratio  of  a  typical  wavelength  relative  to  the  length  scale  over 
which  a  wave  is  modulated.  We  assume  that  dimensionless  variables  can  be 
chosen  in  the  original  problem  leading  to  (2.1)  so  that  both  the  wave 
amplitude  and  the  wavelength  are  0(e). 

Mow  we  describe  how  to  calculate  the  functions  in  (2.2).  The  mean 
solution  u  satisfies  (2.1)  to  0(e).  For  u  regular  in  e  we  expand  it  as 

(2.3)  u  -  u(0,(x)  +  ev(x)  +  0(e2)  . 

We  use  (2.3)  in  (2.1),  expand  A^  and  B  in  Taylor  series  about  u  =  u^°^. 


and  formally  equate  to  zero  the  coefficients  of  and  e1.  It  follows 

that  u<°>  is  a  solution  of  (2.1),  and  v(x)  satisfies  the  system  of 
equations  that  results  from  linearizing  (2.1)  about  u  »  u^. 


(2.4,  l  A(i)(x,u(0),;  *  l  A(i>(x,u(°W0)  +B  <x,u<°^-0  . 

i»1  Xi  i=1  u  Xi  u 


The  notation  A  v  means  £  A  v  .  we  suppose  u^(x)  is  a 

j»1  Uj  3 

function,  and  then  v(x)  is  found  from  (2.4). 


known 


The  sum  of  terms  a^R^*  in  (2.2)  is  a  superposition  of  high  frequency 
waves.  The  jth  wave  depends  on  the  rapidly  varying  phase  4>^3  * (x)/e,  and  the 
wave  is  proportional  to  the  vector  R^(x)  determined  by  the  phase  function 
The  scalar  a^ 3  ^  (x,$ ^ 3  *  (x)/e )  is  the  amplitude  of  the  wave  (if  R^ 3  ^ 
is  normalized  to  unit  length ) . 

The  mean  of  the  wave  amplitude  over  the  rapidly  varying  phase  is  required 


to  be  zero.  That  is 


(2.5) 


Urn  a(j)(x,0)d6)  -  0  . 


In  (2.5)  the  mean  is  taken  over  (0 ,°°) ,  but  it  could  equally  well  be  taken 
over  any  other  semi-infinite  or  infinite  interval. 

The  functions  $^3^,  R^3^  and  a^3^  are  determined  as  follows.  Each 
phase  function  satisfies  an  eiconal  equation 


(2.6) 


det(I  ♦  <i>A<1W°)>)-0  . 


X 

i-1  i 


Notice  that  (2.6)  really  is  an  equation  for  because  A^(x,u)  is 

evaluated  at  u  *  u^®^(x). 

Equation  (2.6)  is  a  first  order  partial  differential  equation  for  a 
scalar.  Therefore  it  can  be  written  as  a  system  of  ordinary  differential 
equations  along  a  set  of  curves  in  R°.  These  curves  are  the  rays  or 
bicharacteristics  of  (2.1)  which  correspond  to  the  solution  u  =*  u^°^  and  the 
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phase  function  $  ^  ^ .  We  introduce  ray  coordinates  ( s  . ,  8  .  3 .  1 ) , 

3  3 » 1  3 1 n- 1 

where  Sj(x)  is  a  function  of  arclength  along  a  ray  and  3j(x)  = 

(3.  .,•••, 8.  . )  is  constant  on  each  ray. 

3,1  3,n-1 


(2.7) 


Then  [8]  the  rays  are  solutions  of 
dx 

1 1 

6j 


ds , 


L(j)(x)A(l)(x,u(0))R(j)(x)  (i  -  1, 


,n) 


In  (2.7)  l/^  and  ^  (the  vector  appearing  in  (2.2))  are  the  left  and 


right  null  vectors  of  the  matrix 


(2.8) 


(D,(i> 


I  C'a 

i=1  i 


In  (2.8),  and  in  the  future,  we  omit  showing  ejq>licitly  that  A^^  and  B 
are  evaluated  at  u  =  u^0^.  Since  (2.1)  is  strictly  hyperbolic  l/  ^  and 
r^7)  are  uniquely  determined  up  to  a  scalar  multiple. 


(2.9) 


The  amplitude  of  the  jth  wave  is  given  by 

a(j)  =  F  (0  /«).(.  ,•)  . 


In  (2.9)  F j  is  an  arbitrary  function  with  zero  mean,  which  describes  the 
wave  profile  by  its  dependence  on  C^/e.  The  factor  E j  in  (2.9)  describes 
the  modulation  of  the  wave  amplitude  due  to  inhomogeneities  and  changes  in  the 
ray  geometry.  It  is 


(2.10)  E.  -  exp  {-/  j[  l  L(j)AU)R(  j)u(0>  +  L<j)B  R(j)]ds.}  . 

3  s0  i«1  Xi  Xi  u  3 


The  integral  in  (2.10)  is  taken  with  respect  to  Sj  along  a  ray  8^  = 
constant.  The  integral  is  taken  from  an  arbitrary  point  s^  =*  sQ(8j)  on  the 
ray. 


The  wavefronts  of  (2.9)  are  given  by  ■  constant.  The  modified  phase 
function  C^(x,E)  is  defined  implicitly  by 


(2.11) 


Cj  -  4>(j)  -  eFj(6j,5j/E)Ij  “  EKj  * 
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In  (2.11)  Ij(3j,Bj  and  Kj(Sj,B^)  are  given  by 


i.  =  I  r  L‘ivi,»,i,E«,»«,».di. 

3  i=1  SQ  xi  u  33 


(2.12) 


i=1  s0  *i  “  3 


The  integrals  in  (2.12)  are  taken  along  a  ray  as  in  (2.10). 

The  expression  (2.11)  for  £ ^  involves  only  Fj  and  does  not  depend 
on  F^  for  any  k  ^  j.  The  different  high  frequency  waves  are  completely 
uncoupled. 

If  (2.1)  is  linear  then  £_.  *  and  the  wavefronts  are  given  by 

=  constant.  Equation  (2.11)  expresses  the  fact  that  the  wavefronts 
£_.  *  constant  differ  from  level  surfaces  of  because  (2.1)  is 

nonlinear.  The  term  £F. I.  in  (2.11)  gives  the  correction  to  the  wavefronts 

j  ] 

because  of  the  waves  nonlinear  self-interaction.  The  term  eie^  gives  the 
effect  of  the  departure  of  the  mean  solution  u(x,e)  from  u^03(x). 

Although  both  the  correction  terms  in  (2.11)  are  small  in  magnitude  the 
nonlinear  case  is  qualitatively  different  from  the  linear  one.  In  the 
nonlinear  case  (2.11)  defines  £j(x,e)  implicitly,  and  because  of  this  £^ 
can  become  a  multivalued  function  of  x.  Then  a  discontinuity  or  shock  must 
be  introduced  into  the  solution  for  a^3  in  order  to  make  a^3  a  single¬ 
valued  function  of  x. 

For  a  conservative  system  such  a  shock  satisfies  an  equal  area  rule. 

Suppose  that  the  ray  fL  =  constant  meets  a  shock  in  the  jth  wave  at 

8^  *  f(8^).  Denote  the  two  wavefronts  meeting  the  shock  at  this  point  by 

£,  =  z  (0.,e)  and  £.  *  z  (8.,e),  Then  z.  and  z-  are  related  by  the 
)  *  ]  )  ^  ^ 

equal  area  rule 
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(2.13)  /  F.(S.,C/e)dC 

Z1  J  3 


-  (z2-z1){Fj(6j,z1/e)  +  Fj(0j,z2/e)} 


We  also  have  from  (2.11)  and  the  definition  of  z^  and  z2  that 
(2.14)  zk  =  $(j)  -  eFj(0j,zk/e)Ij  -  eK..  (k  =  1,2)  . 

All  functions  of  x  in  (2.14)  are  evaluated  at  s^  =  f(P,).  Equations  (2.13) 
and  (2.14)  provide  three  equations  for  the  unknown  functions  z.,  z2  and 
f.  These  equations  are  essentially  the  usual  equations  for  weak  one¬ 
dimensional  shocks  as  described  in  [9] . 

Now  we  briefly  describe  how  to  use  these  results  to  solve  an  initial¬ 
boundary  value  problem.  Suppose  we  are  given  boundary  values  for  a  small 
amplitude  high  frequency  perturbation  about  a  solution  u^°^(x)  of  (2.1) 
which  is  assumed  known.  For  example  in  gas  dynamics  u^°*  might  correspond 
to  gas  at  rest  with  constant  density. 

We  choose  dimensionless  variables  in  which  the  perturbation  amplitude  is 
0(e)  and  the  frequency  is  0(1/6).  Then  v(x)  is  found  by  solving  (2.4),  a 
linear  system  of  partial  differential  equations,  subject  to  boundary  values 
obtained  by  taking  the  mean  of  the  perturbation's  boundary  values  over  the 
fast  phase  variables.  The  infinite  or  semi-infinite  interval  over  which  each 
fast  variable  is  averaged  may  be  choosen  so  as  to  make  the  boundary  values  for 
v  as  simple  as  possible. 

Next  the  m  phase  functions  $^(x)  are  found  by  solving  (2.6),  which 
can  be  done  by  integrating  a  system  of  ordinary  differential  equations. 

Initial  conditions  for  are  chosen  so  that  takes  the  same  value 

on  the  boundary  as  the  appropriate  rapidly  varying  function  on  which  the 
boundary  values  depend.  Once  is  known  it  is  an  algebraic  problem  to 

calculate  the  corresponding  null  vector  . 
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The  amplitude  a^  is  given  by  (2.9)  and  the  function  F^  is  chosen  to 
make  (2.2)  agree  with  the  boundary  conditions.  The  function  E  j  is  given  in 
(2.10)  and  the  modified  phase  variable  is  found  from  (2.11).  If  shocks 

form  they  are  fitted  into  the  solution  using  (2.13)  and  (2.14).  This 
completes  the  solution  of  the  boundary  value  problem.  As  an  example  of  this 
procedure,  a  general  initial  value  problem  for  the  unsteady  one  dimensional 
gas  dynamics  equations  is  solved  in  section  6. 

There  are  two  limitations  on  the  use  of  (2.2).  Firstly  the  results  above 
apply  to  the  nonresonant  case.  Resonant  interaction  between  the  jth  wave  and 
the  kth  wave  (j  /k)  can  occur  if  there  are  nonzero  scalars  c^(x)  and 
ck(x)  such  that  the  vector  p(x)  defined  by 

(2.15)  p  -  CjV$^*  +  , 

satisfies  the  eiconal  equation 

^  /  •  \ 

(2.16)  det{  l  p  AV  }  *  0  . 

i=*1 

The  results  above  apply  provided  (2.15)  and  (2.16)  do  not  hold  for  any 
distinct  j  and  k  at  any  x.  Included  in  this  condition  is  the  case 
p  =  0,  when  and  ^  are  parallel.  This  condition  is  a  sufficient 

condition  for  resonance  not  to  occur.  A  necessary  and  sufficient  condition  is 
derived  in  section  7. 

A  second  limitation  on  the  use  of  (2.2)  is  that  in  general  the  expansion 
will  break  down  for  large  times  of  0(1/e).  One  source  of  this  nonuniformity 
is  the  expansion  (2.3)  of  u.  Another  is  the  cumulative  effect  of  lower  order 
(cubic)  nonlinearities  which  we  expect  to  become  significant  after  times  of 
0(1/e).  if  necessary  one  could  extend  the  validity  of  (2.2)  to  larger  times 
by  introducing  further  multiple  scales  X  =  ex. 


When  (2.1)  is  linear  the  equations  summarized  above  reduce  to  those  of 
geometrical  optics.  We  finish  this  section  by  stating  the  major  effects 
introduced  by  the  nonlinearity  of  (2.1).  There  are  three.  Firstly  for  there 
to  be  no  interaction  between  different  waves  the  mean  wave  profile  (2.5)  must 
be  zero.  Secondly  the  wavefronts  z ^  =  constant  are  distorted  by  the 
nonlinearity  according  to  (2.11),  and  this  can  cause  shocks  to  form.  Lastly 
there  is  the  possibility  of  resonant  interactions  between  waves  when  (2.14) 
holds. 
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3.  Derivation  of  the  Equations 

In  this  section  we  derive  the  equations  for  v(x)  and  a^^  stated  in 
section  2.  Following  Choquet-Bruhat  [1]  we  shall  seek  small  amplitude,  high 
frequency  solutions  to  (2.1)  of  the  form 


(3.1) 


u  *  u  (x)  +  ev(x,0,e) 


In  (3.1)  u'0}  is  a  solution  of  (2.1),  as  before,  and  8  *  (8  ,...,8  )  where 

I  m 

0^  *  $^(x)/e.  The  v  and  are  to  be  determined.  The  difference 

between  (3.1)  and  the  form  of  v  used  by  Choquet-Bruhat  is  that  in  (3.1)  v 
depends  on  m  fast  variables  8^  instead  of  just  one.  This  allows  the 
presence  of  many  waves,  permitting  us  to  deal  with  the  interactions  between 
them. 

We  shall  use  the  method  of  multiple  scales,  in  which  x  and  8  are 
treated  as  independent  variables,  to  obtain  an  asymptotic  expansion  for  v  as 


e  ♦  0.  We  shall  assume  that  v,  v„  and  vfl  are  bounded  functions  of  x 

xi  ei 

and  6  and  that  v(x,8,e)  has  an  asymptotic  expansion  as  t  *  0  of  the  form 


(3.2) 


v(x,8»e)  -  v*°*{x,0)  +  ev^(x,8)  +  0(e2)  . 


We  now  use  (3.2)  in  (3.1)  and  substitute  the  result  into  (2.1).  Then  we 

expand  A^(x,u)  and  B(x,u)  in  Taylor  series  in  powers  of  e  about 

u  »  u^0^,  and  replace  the  partial  derivative  -5 —  by 

xi 

_L  +  1  T  JL 

xi  j-1  xi  j 

Equating  to  zero  the  coefficients  of  e°  and  e1  in  the  resulting  formula 
gives  the  following  equations: 


(3.3) 


l  l  - 


j=»1  i=1  i 


r 


in  n 

»(j)A(i)v(,) 
x  A  v8 

j-1  i-1  Xi  j 


0.4)  i  i 


l  (A(i)"<0)  -  .(i)-(0)-(0) 


i-1 


m  n 


+  A'-'v'-'u  '' ) 
Xi  U  Xi 


♦  I  l  ^'a'^’v'^Bv'01 

“  i.  TV  11  H  11 


j-1  i=1  xi  u 


j 


In  (3.3)  and  (3.4)  A^,  A*^  and  B  are  all  evaluated  at  u  =  u*0^. 

u 

For  a  single  wave  propagating  along  the  j-th  characteristic  only  one  term 
in  the  sum  over  j  in  (3.3)  appears,  and  therefore  it  is  zero.  Motivated  by 
the  idea  of  extending  the  superposition  of  single  waves  from  the  linear  to  the 
weakly  nonlinear  case,  we  shall  look  for  solutions  such  that  each  term  in  the 
sum  over  j  in  (3.3)  vanishes  separately.  Therefore  we  suppose  that 


(3.5)  l  =0  for  j  -  1 . m  . 

i-1  xi  j 

In  order  for  vj0^  to  be  non-zero  the  matrix  in  (3.5)  must  be 
j 

singular.  This  gives  the  eiconal  equation  (2.6)  for  A^(x).  Also  v*0^ 

(j)  3 

must  be  parallel  to  the  right  null  vector  R  J  of  (2.8).  Therefore  we  seek 
a  solution  for  v  ^ 0  ^  of  the  form 


(3.6)  v(0)  -  v(x)  +  l  a(J>(x,0  )R(j)(x)  . 

j-1  3 

Clearly  (3.6)  satisfies  (3.5)  for  an  arbitrary  vector  function  v(x)  and 
arbitrary  scalar  functions  a^^(x,0^).  It  is  not  the  most  general  form  of 
the  solution  to  (3.5),  but  we  shall  see  that  a  solution  of  this  form  can  be 
found  in  the  absence  of  resonance. 

Without  any  loss  of  generality  we  take  v(x)  in  (3.6)  equal  to  the  mean 
of  v^^(x,0)  with  respect  to  0.  The  mean  of  v^®^  is  defined  because 
v^°)  is  a  bounded  function  of  0  at  each  x.  Then  the  mean  of  each 
a^tx,©^)  with  respect  to  0^  is  zero.  Thus: 
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(3.7) 


v(x)  =*  lim  {~  /*•••/„  v^°*(x,9)d9  •••d0  } 

T+m  T*  0  0  1m 

lim  {t  H  a< j)(x,0  )d0  }  *  0  . 


The  functions  v  and  a  ^  will  be  found  from  solvability  conditions  for 

(3.4) . 

The  average  of  the  derivative  of  a  bounded  function  is  zero.  If  the 

expansion  (3.2)  is  uniformly  valid  for  all  x  and  9,  v^°^  and  v^1^  must 

be  bounded  functions  of  x  and  9.  we  use  (3.6)  in  (3.4)  and  average  the 

resulting  equation  over  9  ,...,0  .  The  average  of  the  left  hand  side  is  zero 

1  in 

and  we  obtain  equation  (2.4)  for  the  mean  v(x). 

We  obtain  equations  for  the  m  scalars  a^(x,0_.)  by  using  (3.6)  in 

(3.4) ,  averaging  the  result  over  all  0^  for  k  /  j,  and  then  taking  the 

scalar  product  with  the  left  null  vector  of  (2.8).  The  left  hand  side 

of  (3.4)  is  annihilated  and  we  obtain  an  equation  for  a^.  In  this  equation 
we  recognize  the  sum  of  partial  derivatives  with  respect  to  the  xi  as  a 

directional  derivative  along  the  jth  set  of  rays  (2.7).  we  rewrite  the  sum  of 

these  derivatives  as  a  derivative  with  respect  to  Sj  and  find  that 


satisfies 


(3.8) 


♦  l  {L^V^A^R^R^a^  ♦  L(3Vj)A(i);R<j)}a<l> 

o  «  »  n  y  n  n 


'j  i=1 


xi  u 


xi  u 


+  j/  j  3  )  =  0  . 


In  (3.8) 


(3.9) 


Q™  -  l  (^V^A'^'u'0',  ♦  B  R(^ 
i=1  Xi  U  Xi  U 


The  coefficient  of  a^  in  (3.8)  is  the  propagation  speed  of  disturbances  in 
the  (8^,0^)  plane. 
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The  equations  for  the  scalars  a*^  are  uncoupled.  No  with 

k  ?  j  appears  in  (3.8).  Equation  (3.8)  is  a  nonlinear  equation  for  a(^>  j.n 
two  independent  variables.  It  reduces  to  the  transport  equation  of  (linear) 
geometrical  optics  if  equation  (2.1)  is  linear,  when  the  coefficient  of  aQ^ 

j 

in  (3.8)  is  zero. 

If  the  resonance  condition  (2.4)  holds  then  there  is  a  stronger 
solvability  condition  than  (3.8)  for  (3.4)  and  it  is  not  always  possible  to 
find  a  solution  of  the  form  (3.6)  for  v(0).  This  is  discussed  further  in 
section  7.  Here  we  assume  that  there  is  no  resonance. 

Then  a^  satisfies  (3.8)  which  we  can  integrate  by  introducing 
characteristics.  Let 


(3.10) 


WN'  -  1 

3  3  3  i=1 


vw  *  i  . 

J  1  )  i_,  *i  u 


Then  (3.8)  is 
(3.11) 


VvV "  L(j)2(j)  • 


a^1  +  (M  a(^  +  N.)ai^  +  P  a<j) 
j  3  J 


j 


We  may  write  (3.11)  as 
(3.12) 


j 


oSj  3 


on  the  characteristic  curves 


(3.13) 


—  0 
dSj  j 


M^a 


(j) 


+  N 


The  solution  to  (3.12)  is 

(3.14)  a(j)  -  . 

In  (3.14)  F j  is  an  arbitrary  function  depending  on  6^  and  a  scalar  £ 
which  parametrizes  the  set  of  characteristic  curves  (3.13).  The  function 


Ej  is  defined  by  (2.10) 


We  use  (3.14)  in  (3.13)  and  integrate  keeping  £  constant.  The  result 


is 

S  .  8. 

9j  "  W5)  ^ 0 3  iy«»*j)*j(s»*j)ds  +  V 

(3.15) 

+  5j(0j»C)  * 

In  (3.15)  is  an  arbitrary  function  which  determines  how  £  parametrizes 

the  curves  (3.13).  Usually  we  take  s^(5j«C)  =  £. 

The  final  step  in  the  method  of  multiple  scales  is  to  let 
6^  -  $^(x)/£.  We  use  this  in  (3.15)  and  also  let  £  -  C^/e  in  (3.14)  and 

(3.15) .  Then  we  obtain  (2.9)  from  (3.14)  and  (2.11)  from  (3.15)  after 
multiplying  through  by  6. 

The  only  equation  that  we  still  have  to  derive  is  the  equal  area  rule 
(2.13)  and  we  shall  do  this  in  the  next  section. 


j 

j 


4.  Derivation  of  the  Shock  Conditions 

In  this  section  we  suppose  that  (2.1)  comes  from  a  conservation  law  and 
derive  the  equal  area  rule  (2.13)  for  a^  from  the  Rankine-Hugoniot  shock 
conditions  for  (2.1). 

Therefore  suppose  that  (2.1)  can  be  written  in  conservation  form 

(4.1)  l  {F(i)(x,u)}  +  H(x,u)  -  0  . 

i-1  i 

In  f  .1)  F*^  and  H  are  N-vectors  such  that 


(4.2) 


F^(x,u)  *  A^(x,u)  , 

r  (i) 

H(x,u)  +  2.  F  (*#«)  *  B(x,u) 

i-1  Xi 


The  F(i*  correspond  to  conserved  quantities,  and  H  is  the  source  density 
of  these  quantities. 

The  generalized  Rakine-Hugoniot  shock  conditions  for  (4.1)  which  hold 
across  a  shock  front,  are 


In  (4.3)  [F*1*]  denotes  the  change  in  F*1*  across  the  shock  front,  and 

n  is  the  normal  to  the  shock. 

We  consider  the  jth  wave  with  amplitude  a^* .  Linear  theory  predicts 
that  a  shock  in  this  wave  travels  along  the  linear  characteristics  $^(x)  - 
constant.  This  should  be  the  zeroth  order  result  for  weakly  nonlinear 
waves.  We  suppose  that  the  position  of  a  shock  in  the  j-th  wave  is  given  by 
(4.4)  h(x,0?e)  =  0  . 

From  the  linear  theory  h(x,8>0)  should  be  a  function  of  x  and  0^  alone, 
since  then  the  shock  position  is 
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$^(x)  »  constant  +  0(e) 


Therefore  we  seek  a  power  series  expansion  of  (4.4)  in  the  form 
(4.5)  h(x,0;e)  -  h(0)(x#e_.)  +  eh(1)(x,0)  +  0(e2)  . 

We  assume  that  h^,  h^  and  their  derivatives  are  bounded  functions  of 
x  and  0 . 

The  normal  n  to  the  shock  front  is  proportional  to  Vh  and  its 
component  n^  is  proportional  to 


(4.6) 


h  »^(j)h<0)  +h<°>  +  l  ♦(k,h‘1)  ♦  0(e)  . 

xi  xi  j  xi  k«1  xi  k 


We  use  (3.1)  and  expand  [F^]  in  a  Taylor  series  in  powers  of  e  about 
u  *  u(0).  We  asstime  is  continuous  across  the  shock.  Using  (4.2)  this 


gives 

(4.7)  IFU']  -  SA{i)(x,u(0))[v]  +  j  eVi)(x,u(0))  +  0(e3)  . 

We  us.s  (4.6)  and  (4.7)  in  (4.3)  and  equate  the  coefficients  of  e°  and 
e1  to  zero.  This  gives 


(4.8) 


(4.9) 


.(0)  V  JjJ.U), 

he  l  ♦  ’ 'A  [v  ]  =  0  , 

j  i=1  i 

j  i*i  xi  j  i-i  i 


+  l  h(0)A(i)[v(0)] 
i=1  Xi 


!  !  h^v^v0*])  . 

i=1  k=1  k  i  J 


Equation  (4.8)  is  satisfied  if  [v^0^]  is  proportional  to  Rv-,/.  Therefore 
we  suppose  that  v  and  a^(Mj)  are  continuous  across  the  shock.  Then 


<j>. 


Now  we  use  (3.6)  and  (4.10)  in  (4.9).  We  average  the  resulting  equation 


over  all  8^,  k  ft  j,  and  take  the  scalar  product  with  l/  ^ .  The  result  is 

(4.11)  l  ( [a<^L<j)A(1Vj)h(0)  ♦  1  ((a(^)2]^Vj)A(i)R(jV^ 
i=1  V  Xi  2  xi  u 


In  (4.11)  the  sum  of  partial  derivatives  with  respect  to  can  be  written, 

as  before,  as  a  derivative  with  respect  to  s^.  Then  we  divide  (4.11)  through 
by  [a^-^]  to  obtain: 

(4.12)  h 


^)A(i)R(^R(^)<a  (j)> 
u 


<°>  t  l  />*« 

j  i-1  V.  i 

.h^W^^R^?)  h<°>) 
2  ^  u  0  J 


*  0 


In  (4.12)  we  have  written  [a^2]/2[a^*]  as  <a^*>  where 

(4.13)  <•<»>  («<*>  +  a^}  . 

In  (4.13)  a+^  and  a^^  are  the  limiting  values  of  a^  on  the  two  sides 

of  the  shock,  and  <a^>  is  their  average. 

We  claim  that  for  systems  derived  from  a  conservation  law,  if  S  and 
T  are  any  N- vectors. 


(4.14) 

Proof.  From  (4.2) 


A(1)ST  =  A(1)TS  . 
u  u 


3F(i) 

.(1)  _  1 

*  '  3\  * 


Therefore 
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Comparing  equation  (4.15)  with  (3.8),  we  see  (4.15)  expresses  the  well  known 


fact  that  the  speed  of  a  weak  shock  is  the  average  of  the  propagation  speeds 
on  the  two  sides  of  the  shock.  The  shock  problem  (3.11)  and  (4.15)  always 
involves  only  two  independent  variables,  Sj  and  8^,  however  many 
independent  variables  there  are  in  (2.1),  and  the  equal  area  rule  (2.13)  then 
follows  from  (4.15)  exactly  as  described  in  [9]. 
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5.  Comparison  with  the  Nonlinear lgat ion  Technique 


Laudau  [2]  and  Whitham  [3]  independently  proposed  a  technique  for 
nonlinearizing  geometrical  optics*  In  the  technique,  the  deviation  of  the 
rays  from  their  linear  position  is  neglected,  but  in  treating  the  propagation 
of  disturbances  along  the  rays  the  first  order  nonlinear  correction  to  the 
propagation  speed  is  taken  into  account. 

We  consider  a  single  wave  propagating  on  one  family  of  characteristics 
with  amplitude  a.  Whitham  describes  the  nonlinearization  technique,  as 
follows  [9] . 

Geometrical  optics  provides  a  ray  geometry  and  gives  along  each  ray 

(5.1)  a  »  f(4)E(s)  . 

In  (5.1),  s  is  arclength  measured  along  the  ray.  The  function  E(s)  gives 
the  wave  amplitude,  f  (4)  gives  the  wave  profile,  and  4(b)  is  the  phase 
defined  by 

(5.2)  ‘  -  /;  • 

In  (5.2),  cQ ( s )  is  the  velocity  at  which  disturbances  are  propagated  along 
the  ray  according  to  linear  theory. 

The  nonlinear  velocity  on  the  ray  c(s,a)  is  expanded  for  small  o  as 

(5.3)  c(s,a)  -  Oq(s)  +  c^sja  +  0(o2)  . 

Then  the  result  of  nonlinearizing  (5.1)  is 

(5.4)  a  -  f(T)E(s)  . 

In  (5.4),  the  nonlinearized  characteristic  variable  t  is  defined  implicitly 
by 

<5,5>  ^frr " T  - f(T) /;  ci(8')co2(8,)E(8,)d8’  * 

We  show  that  (5.4)  and  (5.5)  follow  from  weakly  nonlinear  geometrical 
optics.  Let  the  wave  amplitude  be  ea.  Then  a  satisfies  equation  (3.11). 


That  is 


(5.6)  a  +  {M(s)a  +  N(s)}a„  +  P(s)a  =  0  . 

s  o 

In  (5.6),  M(s),  N(s)  and  P(s)  are  given  by  (3.10).  Also  8  =  $/e,  where 
$  is  a  solution  to  the  eiconal  equation  (2.6).  We  may  take  $  to  be  given 
by  (5.2),  provided  the  matrices  are  independent  of  t. 

Suppose  we  linearize  (2.1)  about  u  =  u^0^  and  use  linear  geometrical 

optics.  The  transport  equation  is  obtained  by  dropping  the  terms  in  (5.6) 
proportional  to  the  square  of  the  wave  amplitude.  That  is  we  drop  M( s ) aa0 , 
and  N(s)a0  since  N(s)  is  proportional  to  the  mean  v(x). 

Thus  the  equation  for  a  obtained  by  linear  theory  is 

(5.7)  as  +  P(s)a  =  0  . 

The  solution  to  (5.7)  is 

(5.8)  a  =  F( 8 )E( s)  . 

In  (5.8),  f(8)  is  an  arbitrary  function  and  E(s)  is  given  by  (2.10). 

Let  us  compare  (5.8)  with  the  solution  according  to  the  weakly  nonlinear 
theory.  We  write  (5.6)  as 


(5.9) 


da 

—  +  P(s)a  =  0 
as 


d8 

—  *  M(s)a  +  N( s ) 
as 


on 

(5.10) 

The  solution  to  (5.9)  is 

(5.11)  a  =  F(C)E(s)  . 

We  use  (5.11)  in  (5.10)  and  integrate  taking  our  constant  of  integration  to  be 
£.  Then  we  let  8  =  $/e,  use  (5.2)  and  multiply  through  by  e.  in  this  way 
we  find  that 

(5.12)  t  -  JJ  — -  e?  +  EF(£)  JJ  Hls'lBfsMds' 


+  e  N( s ' )ds' 
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We  see  from  the  weakly  nonlinear  solution  (5.11)  and  (5.12)  that  the 


amplitude  a  has  the  same  form  as  in  tne  linear  theory,  but  the 
characteristic  variable  £  is  constant  along  characteristics  (5.12) 
determined  using  nonlinear  theory.  This  is  just  what  is  supposed  in  the 
nonlinearization  technique. 

When  N(s)  =  0,  (5.4)  and  (5.5)  agree  exactly  with  (5.11)  and  (5.12)  if 

we  let  T  =  e?  and  f(T)  *  £F(T/e),  provided 

(5.13)  M(s)  =  -c1(s)c02(s) 

We  show  that  (5.13)  is  the  case. 

Suppose  u  satisfies 

n  ... 

(5.14)  u  +  £  A  1  (x(u)u  +  B(x,u)  =  0 

i=1  Xi 


The  eiconal  equation  is 

^  /  •  \ 

(5.15)  det{®  I  +  l  *  A  (x,u)}  »  0  . 

C  i«1  Xi 

The  phase  speed  c(x,u)  equals  -4^/ |V®|,  and  from  (5.15)  satisfies 

n  ... 

(5.16)  det{c(u)I  -  l  n.AU,(u)}  =  0  . 

i=1  1 


The  speed  c  also  depends  on  the  direction  of  the  normal  to  the 


characteristic  surface,  n  = 
and  A^^  on  x  in  (5.16). 


V® 

|V*| 


.  We  omit  showing  the  dependence  of  c 


Let  L(u)  and  R(u)  be  the  left  and  right  null  vectors  of 


v  ( i ) 

c(u)I  -  2  (u)*  Then 

i=1  1  n 

(5.17) 


£  n.A(i)  (u)R(u)  =*  c(u)R(u) 


i=1 


We  expand  about  u  =  0  (there  is  no  loss  in  generality  in  taking  u^0^  =  0): 


-22- 


R(u)  =  R  +•  RyU  +  0(u  )  , 

(5.18)  2 

L(u)  =  L  +  L^u  +  0{u  )  , 

X  n.A(i)(u)  =  X  n.A(i)  +  X  T».A(i)u  +  0(u2)  . 

i=1  1  i=1  1  i=1 

On  the  right  hand  side  R,  L,  A* i ^  and  their  derivatives  are  evaluated  at 
u  =  0. 

We  use  (5.18)  in  (5.17)  and  equate  coefficients  of  powers  of  u.  Then 


(5.19) 


l  niA(i)p  =  cqr  , 


i=1 


(5e20) 


(  I  n.A(i)u  -  c  ui)r  =  -(  X  n. 
i=1  1  U  U  i*1 


A(l)  -  cnl )r  u  . 
0  1  u 


(i) 

From  (5.19),  R  is  the  right  null  vector  of  \  n  a'  -  cl. 

i=1 

Similarly  L  is  the  left  null  vector.  Then  we  take  the  scalar  produce  to 
(5.20)  with  L  and  solve  the  resulting  equation  for  cuu.  We  find 

n 

(i) 


(5.21) 


c  u  = 
u 


l  n.LA1  'uR 
i-1  1  U 


LR 


To  find  the  velocity  along  the  linear  rays,  we  fane  h.  =  $  /l^<H  in 

i 

(5.21)  where  $  is  given  by  (5.2).  Then  since  cQ  =  -9t/l^9|  and  =  1 . 

(5.22) 


n.  =  -c  9 
1  ox, 


Next  the  derivative  with  respect  to  s  in  (5.6)  is  taken  along  the  rays 
(5.23) 


~  =  LR  ,  f  x.  =  LA(i)R 
ds  ds  l 


We  differentiate  (5.2)  with  respect  to  s,  and  use  4^  =  0,  since  9  is 

ds 


constant  on  the  rays.  The  result  is 
(5.24) 


dt  -1 
ds  °0 
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Therefore  from  (5.23)  and  (5.24) 

(5.25)  LR  =  c"1  . 

(This  is  the  normalization  of  L  and  R  required  so  that  the  parameter  s 
is  arclength.) 

We  use  (5.22)  and  (5.25)  in  (5.2).  Then 

(5.26)  c  u  =  -cl  l  <p  IA(±)uR  . 

u  0  .  „  x.  u 

1=1  i 

Finally  we  let  u  =  or  in  (5.26),  and  use  the  result  in  (5.18).  Then 

(5.27)  c(u)  =  cQ  -  ac^M(s)  +  0(a2) 

n 

In  (5.27),  we  have  used  (3.10)  to  write  £  $  lAtlJRR  =  M(s).  Now,  c  was 

i=1  Xi  U 

expanded  in  the  nonlinearization  technique  as  (5.3).  We  see  that  (5.13) 
follows  from  (5.3)  and  (5.27).  This  shows  that  for  weak  waves  the 
nonlinearization  technique  follows  from  weakly  nonlinear  geometrical  optics, 
when  we  apply  it  to  single  waves  (N  =  0). 
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6.  Comparison  with  the  Method  of  Characteristics 

Systems  of  hyperbolic  equations  in  two  independent  and  two  dependent 
variables  can  be  solved  by  introducing  characteristic  coordinates  [10].  In 
this  section  we  show  that  the  results  calculated  using  characteristic 
coordinates  for  the  one  dimensional  gas  dynamics  equations  agree  to  0(£) 
with  those  calculated  using  weakly  nonlinear  geometrical  acoustics. 

The  gas  dynamics  equations  for  plane  (N  =  0)  axisymmetric  (N  =  1), 
and  spherically  symmetric  (N  =  2)  flows  are 


(6.1) 


c  +  uc  +  c(u  +  N  — )  =  0  , 

t  x  2  x  x 


u  +  uu  +  — — -  cc  =  0  . 

t  x  Y-1  x 


In  (6.1)  u  is  the  gas  velocity  and  c  is  the  local  sound  speed.  To  be 
definite  we  shall  consider  an  initial  value  problem  for  (6.1),  and  take  as 
initial  values 


c(x,t=0)  -  c0  +  ecQg(x,x/e)  +  0(e2) 

(6.2) 

u(x,t=0)  =  ecQf(x,x/e)  +  0(£2)  . 

In  (6.2)  f(x,6)  and  g(x,0)  are  arbitrary  continuously  differentiable, 

bounded  functions. 

Here  we  shall  solve  (6.1)  and  (6.2)  for  plane  and  spherical  flows.  The 
results  for  axisymmetric  flows  are  entirely  analogous  to  those  for  spherical 
flows,  although  the  formulae  are  more  complicated.  The  details  are  given  in 
(11]  . 

We  shall  solve  (6.1)  and  (6.2),  first  by  weakly  nonlinear  geometrical 
acoustics  and  second  using  characteristic  coordinates.  Then  we  shall  show 
that  the  two  solutions  agree. 
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To  solve  (6.1)  and  (6.2)  by  weakly  nonlinear  geometrical  acoustics  we 


seek  a  solution  of  the  form  (A. 16)  with  Cg(x,t)  =  Cq  and  U(x,t)  =  0.  In 
this  one-dimensional  problem  there  are  only  sound  waves.  The  phase  functions 
4>^(x,t)  satisfy  (A. 7)  which  is 

(6.3) 


(<j))2-cXj,.2=0 


We  want  the  fast  variables  4>^*/e  in  (A.  16)  to  equal  the  fast  variable  x/e 


appearing  in  (6.2)  at  t  =■  0.  Therefore 

(6.4) 


(J>(j)(x,t=0) 


There  are  two  solutions,  to  (6.3)  and  (6.4): 


$(±>(x,t)  -  x  +  CQt 


(6.5) 

From  (A. 18)  the  corresponding  null  vectors  are 


(6.6) 


K(1)  -  C. 


■  ’  1 
±2/(Y-1)J 


where  we  have  chosen  the  arbitrary  scalar  X±  in  (A. 18)  to  be 
(6.7)  X±  *  *2c0/(Y-1)  . 

Using  (6.5)  and  (6.6)  in  (A. 16)  we  have 


(6.8) 


•  - 

r  i 

c 

c0 

c(x,t) 

# . .  x-c_t 

i 

* 

+  e 

+  tea  (x,t,  ) 

u 

0 

u( x, t ) 

u  t 

2/(Y-1) 

—  4 

•»  _ 

r  ^ 

-  -4 

,  ,  x+c„t 
+  ec0a(  (x,t, — - — ) 


1 


■2/ ( Y~1 ) 


+  0 (e2) 


The  means  c  and  u  satisfy  (A. 17)  with  U  =  0 

Y-1 

ct  " 

(6.9) 


c  +  c.(u  +  N  — ) 
t  2  0  x  x 


u.  +  — —  c  c  =  0 
t  Y-1  o  x 

Initial  conditions  for  (6.9)  are  found  by  averaging  the  initial  conditions 
(6.2)  for  (6.1): 
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(6.10) 


c(x,t=0)  =  g(x)  , 
u(x,t=0)  =  f(x) 


In  (6.10)  we  have  defined 

(6.11)  f(x)  =  lim  tj  /T  f(x,9)d0>,  5(x)  =  lim  /J  g(x,0)d&}  . 

T+40  x+0 

For  N  =  0  the  solution  to  (6.9)  and  (6.10)  is 


(6.12) 


o  =  p0(x-e0t)  +  q0(x+c0t)  , 

®  ■  7TT  {P0(x“cot)  "  Vx+cot,}  ' 


where 


(6.13) 


For 


-  CQ  V— 1  _  _ 

P0(x>  *  2  ^  2~  f<x)  +  9(x)>  » 

qQ(x)  *  {-(~0f(x)  +  g(x) }  . 

N  *  2  the  solution  to  (6.9)  and  (6.10)  is  [9] 


P2(*”cot)  +  ^2<X+C0t) 


(6*14) 

2  f*2(x"C0t)  "  «2{*+C0t}  *2U~C0t}  ‘  Q2(x+C0t), 

U  y_ i  I  x  *  2  *  ' 

x 

where 


P2(x) 

(6.15) 


d_ 

dx 


P2(x) 


*2(x)  *  d^52(x) 

Now  we  calculate 
4»(t)  are 


xf(x)  +  /X  ffx'Jdx1  +  xg(x)}  , 

{-(^Jxffx)  -  (^-)  /Xf(x')dx*  +  xg( x ) }  . 

Ray  coordinates  (s±,8^)  corresponding  to 


(6.16) 


s.  3  x  ±  c.t 


0^  -  x  +  c„t 


With  ray  coordinates  (6.16)  we  must  take  in  (A. 8)  to  be 


(6.17) 


From  (6.5)  we  have 


"j'-°o2/2  • 


(6.18) 


,(±)  2..(±)  VT  2 

0  “  c  A®  —  -Nc 

ytt  0  Y  0‘ 


Ncq/x  . 


We  use  (6.17)  and  (6.18)  in  (A. 10)  then  integrate  with  respect  to 


s^  =  x  ±  cQt  along  the  ray  3^  =  x  +  cQt  =  constant.  This  gives 


(6.19) 


E(t>  .  „-»/2 


Equation  (6.19)  is  the  usual  result  of  linear  geometrical  acoustics:  the 

u 

cross-sectional  area  of  a  ray  tube  increases  like  x  ,  so  the  wave  amplitude 
decays  like  x-N//^  in  order  that  energy  (which  is  proportional  to  the 
amplitude  squared)  is  conserved  along  the  ray  tube. 


We  use  (6.19)  in  (2.9)  and  obtain 


(6.20) 


a(±)  =  f^N)(x  +  cQt,  C±/e)x 


The  functions  F^  are  found  from  the  initial  conditions.  The  fast  variable 
C±/e  in  (6.20)  should  equal  the  fast  variable  x/e  in  (6.2)  at  t  =  0. 


Therefore 


(6.21) 


C±(x,t-0,e)  »  x 


Then  we  use  (6.20)  in  (6.8),  substitute  the  result  into  (6.2),  use  (6.21),  and 
solve  for  F^1^.  This  gives 


f^N> (x,?/e )  =  xN/2  Y  {f (x,?/e)  -  f(x)}  +  — -  {g(x,?/e)  -  g(x)}  , 


(6.22) 


F_  ;(x,G/e)  =  -xw/‘  Y  {f (x,5/e)  -  f(x)>  +  ~~  {g(x,C/e)  -  g(x)}  . 


We  use  (6.5),  (6.7),  (6.16),  (6.17)  and  (6.19)  in  (A. 11)  and  integrate. 


To  satisfy  (6.21)  we  take  8^(3^)  “  8^  80  that  =*  t  =  x  at  t  *  0, 

when  s^_  =  3^.  Thus  we  find 


(6.23) 


,(±) 


<*  ■  - c  ; it;  4.=u0ti““|dsi  • 


(N), 


„  X+C_t 
E  C  0 


In  (6.23) 


(6.24) 


C-*c„t  . 


h<2) .  10,  _* 


x+c0t 


The  integral  of  c  and  u  in  (6.23)  is  taken  with  respect  to  s±  keeping 
6^  constant. 

This  completes  the  solution  of  (6.1)  and  (6.2)  by  the  method  of  weakly 
nonlinear  geometrical  acoustics.  The  solution  is  given  by  (6.8).  In  (6.8) 
c  and  u  are  given  by  (6.12)  and  (6.13)  (N=0)  or  (6.14)  and  (6.15)  (N=2). 

The  amplitudes  a*^  are  given  by  (6.20)  and  (6.22)  -  (6.24). 

Now  we  solve  the  same  problem  using  characteristic  coordinates.  When 
written  with  respect  to  characteristic  coordinates  A  and  B  (6.1)  become 
(8] 


2 

- 7  C  -  u  + 

Nuc 

o 

II 

4J 

,  X,  = 

(u-c)t  , 

Y— 1  A  A 

X 

A 

A 

A 

c  +  u  + 
Y-1  B  B 

Nuc 

X 

tB=  ° 

'  *B  = 

( u+c ) t 

B 

(6.25) 


Without  any  loss  of  generality  we  can  impose  that  A=B=x  at  t  =  0.  Then 
from  (6.2)  the  initial  conditions  for  (6.25)  are 


c(A,B=A)  =  Cq  +  Ec0g(A,A/E)  +  0(E  )  ,  x(A,B=A) 


(6.26) 


u( A , B=A)  =  ecQf (A,A/e)  +  0(e  ) 


,  t(A,B=A)  =  0  . 


We  shall  solve  (6.25)  and  (6.26)  by  seeking  a  generalized  asymptotic 
expansion  for  c,  u,  x  and  t  of  the  form 
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(6.27) 


( 1 )  2 
c  =  cQ  +  £c  (A,B,£)  +  0 ( e  )  , 

( 1 )  2 
u  =  £u  (A,B,£)  +  0(e  )  , 

x  =  x(0)(A,B,£)  +  ex(1)(a,b,e)  +  0(e2) 

t  =  J0)(A,B,£)  +  £t(1)(A,B,E)  +  0(£2) 


We  use  (6.27)  in  (6.25)  and  (6.26),  and  equate  explicit  powers  of  £.  From 
the  coefficients  of  we  find  that 


(6.28) 


x<°>  +  'n^0)  =  0 
A  0  A 


x<0)  -  JO)  .  0 

B  OB 


X<0) (A,B=A)  «  A 


t(0)(A,B=A)  *  0 


The  solution  to  (6.28)  is 


(6.29) 


K(<,>  -  l  (»«, 


v‘0>  - 


Equating  the  coefficients  of  e  and  using  (6.29)  we  find  that  u 


(1) 


and  J1)  satisfy 


(6.30) 


'  =  0  #  c(  1  )  (A,B=A,£)  -  cQ(g(A,A/E)  , 
+  U^1*  +  N^~b  =  0  ,  u(  1  )  (A,B=A,E)  =  cQf  (A, A/E)  . 


Also  x* 1 )  and  J1)  satisfy 


(6.31) 


+  Vi”  =  {-c<1)  +  u(1,}/2co 


XB1}  -  V*"  =  tc(1>  +-{1)>/2c0 


x^ 1 ^ (A,B=A,E)  =  0  , 

t(1) (A,B=A,E)  =  0  . 


Equations  (6.30)  have  the  same  form  as  the  linearized  gas  dynamics 
equations,  so  we  can  write  down  their  solution.  For  plane  flows  (N*0)  the 
solution  to  (6.30)  is 
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(6.32) 


C(1)(A,B,£)  =  pQ(A,A/£>  +  (B,B/£  )  , 

(  1 )  2 

U  ' (A,B,£ )  =—  {pQ (A,A/E )  -  qQ (B,B/£ ) }  , 


where 


P0<A,A/e)  =  ■—  f(A,A/e)  +  g(A,A/£ ) }  , 


(6.33) 


q0(B,B/e)  =  y  {-(^y)f  (B,B/£)  +  g(B,B/E)}  . 

For  spherically  symmetric  flows  (N=2)  the  solution  to  (6.30)  is 


(6.34) 


c  ' (A,B,£ )  =—  {p2(A,A/e,e)  +  q2(B,B/e,e)} 


(  1 )  2  2 

u  '  (a,b,£ )  =  7TTfA+i(p2(A'A/E'£)  -  q2(B,B/e,e)> 


(A+B) 


2  (P2(A,£)  -  Q2(B,£))}  . 


In  (6.34) 


(6.35) 


P2(A.A/£,£)  =|^P2<A#£)  Af(A,A/£) 


+  *y  /A  f  (A1  »A' /£ )dA’  +  Ag(A,A/£ )  }  , 


q2(B,B/e,e)  =  fg  Q2<B,e)  =  y  {-(^)Bf  (B,B/£ ) 


/B  f  (B'  .B'/EjdB*  +  Bg(B,B/E )  }  . 


We  integrate  (6.31)  to  find  x^ 1  ^  and  t^1^  in  terms  of  c^ 1 ^  and 
u(1>,  and  use  the  result  and  (6.29)  in  (6.27).  After  rearrangement,  this 
gives  the  following  expressions  for  A  and  B 


(6.36) 


A  -  x  -  c.t  -  vr~  /»  {c(1)(A,B',e)  +  u(1)(A,B',e)}dB'  +  0(e2)  , 

□  2Cq  A 

B  =  x  +  C  t  +  ~~  /B  {c(1  }  (A',B,e)  -  u^^(A'»B,£)  }dA'  +  0(£2)  . 

0  2c0  A 


In  (6.36)  and  1 ^  are  given  by  (6.32)  when  N  =  0,  and  by  (6.34) 

when  N  *  2. 

This  completes  the  solution  of  (6.1)  and  (6.2)  using  characteristic 
coordinates.  The  functions  c  and  u  are  given  by  (6.27)  and  (6.32)  or 
(6.34) ,  while  x  and  t  are  given  by  (6.36).  Now  we  show  that  this  solution 
equals  the  solution  found  using  weakly  nonlinear  geometrical  acoustics  to 
first  order  in  e. 

We  shall  need  the  following  lemma. 

Lemma  6.1.  Let  f(x,0)  be  a  differentiable  function  of  x  and  integrable 
with  respect  to  0  on  [0/-).  Define 

(6.37)  iu)  -  iim  4  fu,e>ae>  , 

■p*00  1 

and  suppose  that  f(x)  is  differentiable  and 

(6.38)  f  (x)  =■  lim  4  JT  f  (x,0)d0} 

x  r  u  x 

uniformly  in  x.  Then 

(6.39)  lim  fb  f(x,x/e)dx  =  Jb  f(x)dx  . 

e+0  a  a 

Proofs  Let 

(6.40)  F(x,9)  =  f(x,0)  -  f ( x )  , 

and 

(6.41)  u( x,e )  -  /*/E  F( x, 0 )d0  . 

Then 

(6.42)  ux  “  e  F(x'x/e)  +  / q/C  Fx(x,0)d0  . 
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We  multiply  (6*42)  by  e  and  integrate  with  respect  to  x  over  [a,b] •  This 
gives 

(6.43)  /b  F(x,x/e)dx  =  e  /b/£  F(b,0)d0  -  e  /a/£  F(a,0)d0  - 

a  o  0 

-  e  Jb  dx  f*/C  F  (x,0)d0  . 

a  0  x 

Taking  the  limit  of  (6.43)  as  e  +  0  and  using  (6.37),  (6.38)  and  (6.40),  it 
follows  that 

(6.44)  lim  /b  F(x,x/e)dx  =0 

£♦0  a 

Then  (6.39)  follows  from  (6.40)  and  (6.44).  ■ 


In  fact,  with  certain  additional  assumptions  on  f(x,9)  it  is  shown  in  [11] 
that  as  e  0 


/b  f(x,x/e)dx  =  /b  f ( x ) dx  +  0(e) 
a  a 

We  shall  also  need  the  following  three  equations. 


From  (6.16)  and  (6.36) 


(6.45) 


x  -  cQt  =  s_  =  $+  =  A  +  0(e)  , 

x+cQt=s+=  8_  =B+0(e) 


Then  if  h(x,y,0,a)  is  any  continuously  differentiable  function  of  x  and 
y,  using  (6.45)  and  expanding  in  a  Taylor  series  about  A  =  x  -  cQt, 

B  =  x  +  cQt  we  have 

(6.46)  h(A,B, A/e ,B/e )  =  h(x-c0t,x+cQt,A/e,B/e)  +  0(e) 

Finally  suppose  two  continuously  differentiable  functions  w(A,B)  and 
w(x,t)  are  related  by 

w(A,B)  =  w( x , t )  +  0(e)  . 

Then  using  (6.45)  and  expanding  in  Taylor  series  we  have 
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(6.47) 


x+c„t 


Jj  w(A/B' )dB'  =  /x_c°t  wds+  +  0(e)  , 


,B  -  ,x+c0t  - 

w(A',B)dA'  «  /  wds  +  0(e)  . 

A  x"c0t 


In  (6.47)  the  integral  with  respect  to  s+  is  taken  keeping  0+  constant  and 
the  integral  with  respect  to  s_  is  taken  keeping  0_  constant. 

Now  we  compare  the  two  solutions.  We  take  the  plane  and  spherical  cases 


separately  and  consider  plane  flows  first.  Using  (6.13),  (6.22)  and  (6.33)  we 

can  rewrite  the  solution  (6.32)  for  c^1^  and  u^ 1 ^  as 

c(1)(A,B)  =  c (A, B)  +  cqf|0) (A, A/e)  +  cqf[0) (B,B/e)  , 


(6.48) 


In  (6.48) 


u(1)(A,B)  =  u( A, B)  +  F^°  }  (A, A/E)  -  F^0)(B,B/E)  . 


(6.49) 


c(A,B)  -  pQ(A)  +  qQ(B)  , 


i(A,B)  =  {Pq  (A)  -  qQ  (B)  }  . 


Now  using  (6.12)  and  (6.46)  in  (6.49)  we  have 


(6.50) 


c (A, B)  =  c ( x , t )  +  0(e)  , 

u(  A ,B )  =  u(x,t)  +  0(e) 


We  use  (6.46)  and  (6.50)  in  (6.48)  to  obtain 

c(1)(A,B)  =  c(x,t)  +  cqf|°  ^ (x-cQt,A/e)  +  c0F^0)(x+c0t,B/e)  +  0(e) 

(6‘51)  2c  2c 

u(1>(A,B)  =  u(x,t)  +  — J  f|0> (x-cQt,A/e)  -  ~  F^0) (x+cQt,B/e)  +  0(e)  . 

Using  (6.51)  in  (6.27)  and  comparing  the  result  with  (6.8)  and  (6.20)  we  see 
that  the  two  solutions  for  c  and  u  agree  to  0(e)  provided  that 


(6.52) 


A/e  =  c+/e  +  0(e),  B/e  =  c_/e  +  0(e)  . 


To  show  that  (6.52)  does  hold,  we  use  (6.48)  in  (6.36),  and  integrate. 


This  yields 
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(6.53) 


-  «-c0t  -  «  lii  ♦  f  £1  /f 


-  fB  {c(A,Br )  +  u(  A  ,  B ' ) }dB*  +  0(e2)  , 

2c0  1 A 

B  =  x+cot  +  e  lil  Fl0)(B,B/e)  Sli  -  |  3=1  f*  r;#,(A'l*V«m' 


+  /B  { c( A 1 , B )  -  u(A*,B)}dB’  +  0(e2)  . 

2c0  ;  A 

Now  from  (6.22)  F^N)(x,0)  has  zero  mean  over  6,  and  therefore  using 


lemma  6 . 1 


(6.54) 


fB  f!N) (A’fA'/eJdA*  =  0(E)  , 

A  + 

JB  f(N) (B' ,B,/e)dB’  =  0(£)  . 

'  A  “ 


Also  using  (6.50)  and  (6.47)  with  w  =  c  ±  u  we  have 


x+c  t 
r  0 


r  {c(A,B'  )  +  u( A,B'  )  }dB*  =  /  {^+iI}dsA  +  0(e) 

A  X—C^t  + 


(6.55) 


rx+c0t 


ff  {c(A',B)  -  u(A' ,B) }dA'  =  /  U.  {c-u}ds  +  0(e) 

A  X“Cgt  ” 


We  use  (6.46),  (6.54)  and  (6.55)  in  (6.53).  This  gives 


x+o.t 
e  r  0 


A  =  x-c0t  -  e^F )  (x-c0t,A/e)c0t  -  —  /x_C()t<c+a}ds+  +  0(e  ) 


(6.56) 


x+c  t 

B  =  x+cQt  +  e  f^° ^  (x+cQt , B/e )cQt  +  ~  /x_c°t^-u}ds_  +  0(e2)  . 


Comparing  (6.56)  with  (6.23),  and  using  (6.5)  and  (6.24),  shows  that  (6.52)  is 
true  for  plane  flows.  This  completes  the  confimation  that  the  two  methods 
agree  for  plane  flows. 


Next  consider  the  spherical  case.  Applying  lemma  6.1  to  (6.35)  we  have 


P2(A,A/e,e)  =  p2(a,c)  =  ^  Af (A, A/e )  + 


+  /A  f(A')dA'  +  Ag(A,A/e  ) }  +  0(e)  , 


(6.57) 


q2(B,B/e,e)  =|^Q2(B,e)  =  {-(^)Bf(B,B/e)  - 

-  (-^-)  fB  f(B,)dB'  +  Bg(B, B/e ) }  +  0(e)  . 

Integrating  (6.57)  to  obtain  P2  and  Q2  and  applying  Lemma  6.1  again  gives 

c 

P2(A'E)  “  ~2~  ^  A,f(A’>dA'  +  ~  fA  f (A")dA"  + 

+  /A  A1 g(A' )dA' }  +  0(e)  , 

c 

(6.58)  Q2(A,e)  =  {-(-^O  fB  B 1  f  ( B '  )  dB '  -  (^-)  fB  f(B")dB"  + 

+  /B  B'g(B' )dB'}  +  0(e)  . 


Thus  from  (6.58) 


(6.59) 


P2(A,e)  =  P2(a)  +  0(e)  , 

Q2<A,e)  =  Q2(a)  +  0(e)  , 


where  from  (6.15) 

(6.60) 

Now  using  (6.15)  and  (6.22)  we  can  rewrite  (6.57)  as 


da  Vs'  "  ?2IM'  “  52<b' 


( 2 ) 


(8.61) 


P2(A,A/e,e)  =  p2(a)  +  cqf^  '(a, A/e)  +  0(e) 


-  (  2 ) 

q2 (B,B/e ,e )  =  q2(B)  +  cQF_  ' ( B , B/e )  +  0(e) 


Then  we  use  (6.59)  and  (6.61)  in  (6.34)  which  gives 
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(6.62) 


2c, 


c(1)  =  c(A,B)  +  ~  {f|2)(A,A/£)  +  F^2)(B,B/£)}  +  0(e)  , 


2c 

u(1)  =  u(A,B)  +  ^j-  ^  {f|2)(A,A/C)  -  F^2)(B,B/e)}  +  0(e)  . 


In  (6.62) 


(6.63) 


c(A,B)  =  —  {p2 (A)  +  q2(B)}  , 


*(A'B)  =  FT  15 I  {p2(a)  “  VB)}  "  FT  ^2  {iVA)  "  a2(B» 


We  use  (6.46)  and  (6.14)  in  (6.63)  which  leads  to  (6.50).  Then  we  use 
(6.46)  and  (6.50)  in  (6.62)  and  find 


C(1)  =  c(x,t)  +  {f|2^ (x-cQt,A/e)  +  F^2* (x+cQt,B/e)}  +  0(e)  , 

(6.64) 

u(1)  =  u(x,t)  +  -2—  -2.  {f*2) (x-c  t,A/e)  -  F(2> (x+c  t,B/e)}  +  0(e) 

T  «  X  *  U  —  U 


Using  (6.64)  in  (6.27)  and  comparing  the  result  with  (6.8)  and  (6.20)  we  see 
that  the  two  solutions  agree  provided  (6.52)  holds. 

To  show  that  (6.52)  holds  for  spherical  waves,  we  use  (6.62)  in  (6.36) 
which  gives 


=  x-c0t  - 


.  Y+1  (2) 

e  — -  F 
Y-1  + 


.*.AA)log(|2)  .  c  !l 


3-Y 


(2) 

„  FU,(B',B'A) 
r  B  - 


A+B* 


dB* 


-  /B  {c(A,B')  +  u(A,B')}dB'  +  0(e)  , 

B  -  x.=0t  -  c  lil  .  .  (£1)  /B  -  «. 

+  /?  {c( A ' , B )  -  u(A',B)}dA'  +  0(e)  . 

0 


Now  we  use  (6.46),  (6.54)  and  (6.47)  with  w  =  c  ±  u  in  (6.65).  The 
result  is 
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Y+1 

A  =  x-cQt  -  £ 

F+  (x-c0t,A/e)log(x_c^t) 

e 

“  2co 

x+c0t 

x~cot 

{c+u}ds+  +  0(e)  , 

(6.66) 

Y+1 

B  .  x+c0t  -  e  — 

F-  (x+cOt'B/£)log^+c0t^ 

e 

x+c0t 

x~c0t 

{c-u}ds  +  0(e) 

Comparing  (6.66)  with  (6.23),  and  using  (6.5)  and  (6.24),  shows  that  (6.52) 
holds  for  spherical  flows.  Thus,  weakly  nonlinear  geometrical  acoustics  and 
the  method  of  characteristics  agree  to  first  order  in  the  wave  amplitude  e. 


7.  Resonance  conditions 


Equations  for  v(x)  and  a^(x,0)  were  found  in  section  3  as 
solvability  conditions  for  (3.4).  In  this  section  we  show  that  there  can  be 
solvability  conditions  additional  to  those  considered  in  section  3,  and  that 
the  solution  derived  in  section  3  satisfies  all  these  solvability  conditions 
whenever  the  resonance  condition  (2.15)  does  not  hold  for  distinct  j  and 
k  at  any  x.  Even  if  (2.15)  does  hold,  and  there  is  the  possibility  of 
resonance,  the  solution  derived  in  section  3  may  still  be  valid  and  we  give 
necessary  and  sufficient  conditions  for  this  to  be  the  case. 

We  abbreviate  (3.4)  to 

(7.1)  l  l  <j)A(i)v0  =  f  ( 8)  . 

j=1  i-1  i  j 

We  have  dropped  the  superscript  on  v^ 1 ^  and  omit  showing  any  x-dependence 
explicitly,  because  x  is  effectively  constant  in  (7.1).  We  shall  derive 
conditions  which  f(0)  must  satisfy  if  (7.1)  has  a  bounded  solution  for 
v(0). 

In  fact  all  we  require  on  v(8)  is  that 

(7.2)  v(0)  =  o( | 8 | )  as  | 0 |  ♦  •  . 

Then  using  (7.2) 

ev(1)(<fr(1)/e,...,$(m)/e)  =  o{1)  as  e  +  0  # 

asymptotic  expansion  (3.2)  remains  valid  when  6^  is  evaluated  at 
v(8)  satisfies  (7.2)  then 

lim  <?  fo  ve  d9j}  =  0  * 

any  derivative  of  v(0)  has  zero  mean. 


and  the 


If 


(7.3) 


That  is 


-39- 


Therefore  averaging  (7.1)  with  respect  to  6  we  find  as  in  section  3 
that  f { 9 >  has  zero  mean,  which  gives  (2.4)  for  the  mean  v  of  v^°^(0). 
Next  suppose  that  there  is  a  vector  A  =  (A  ,...,Ab)  such  that 


(7.4) 


where  p  *  (p^ 


I  v* 

j=1  3 


(j) 


,...,p  )  satisfies  (2.15) 


11  (  i  ) 

Denote  by  L  the  left  null  vector  of  £  p  a'  .  We  assume  throughout 

i=1  1 

this  section  that  (2.1)  is  strictly  hyperbolic,  so  L  is  uniquely  defined  (up 


to  a  scalar  factor).  Also  let 
(7.5)  o 


m 

l  A  .  0  .  , 


j-1 


3  3 


and  define  the  hyperplane  S(o)  by 

(7.6)  S(O)  »  {0  e  Rm  s  o  »  constant}  . 

Then  we  average  (7.1)  with  respect  to  0  over  S(0) .  To  do  this 

introduce  coordinates  ( ° , n ,,..., n_  , )  in  where  the  n  depend 

I  m- 1  ] 

linearly  on  0.  Then  using  (7.5),  (7.1)  becomes 

m  n  m-1  m  n  3n.  ...  ... 

<7-7’  *  J-,  Vi! A  v« 7  kJ,  *  J,  *r  »  v  *  f<6>  • 

j-1  i-1  J  i  k-1  j-1  i=1  3  i  * 

Now  we  average  (7.7)  over  S (0),  integrating  with  respect  to  n.,...,n  ., 

1  m— 1 

when  the  mean  of  the  terms  proportional  to  v^  is  zero  from  (7.2).  Then  we 

take  the  scalar  product  of  the  result  with  L,  which  makes  the  coefficient  of 

v  zero.  Therefore 
a 


(7.8) 


L(f ) 


avS 


In  (7.8)  is  t*le  avera9e  f  over  S(0). 

<7-9’  <1*  dv-5v.)  • 


Thus  ( f ) 


avS 


T>oo  T 

is  a  function  of  0. 


40 


In  general  it  is  possible  for  there  to  be  many  vectors  A  such  that 

J  AA^a^  has  a  given  left  null  vector  L.  The  set  of  all  such  A  forms 

..ix. 

1/3  i 

a  vector  subspace  of  tP .  Pick  a  basis  {A^,...,Ay}  of  this  space,  where 


A„  =  (A. _,...,A  ),  and  define  a  and  S(o)  analagously  to  (7.5)  and 

x  Jci  xm  x 


(7.6)  s 


—  I  A .  .8 .  (^  -  1,..*,v)  , 

4  j-1  43  3 


(7.10) 


S(o)  =  {0  e  R  s  (<J  ,...,0^)  =  o  =  constant} 


Then  exactly  as  before  we  obtain  the  solvability  condition  (7.8).  We 
obtain  one  such  condition  for  each  subspace  S  of  the  kind  defined  above. 

Next  we  show  that  the  solution  (3.6)  for  v^  always  satisfies  these 
solvability  conditions  if  (2.15)  never  holds.  Using  (3.6)  in  (3.4)  and 
comparing  with  (7.1)  we  obtain 

n 


(7.11)  -f (x,6 )  «  l  {A(i)vx  +  A^i)vux°>}  +  Buv 
i=1  i  i 


m  n 


+  I  I  {AV1,R 

k=1  i=1 


(i)„(k)  (k)  .  . (k)  (i)~  (k)  (k)  (k)  (k), 

"a  +  <|>AvRaa+Qa} 
x.  x,  u  0, 

i  i  k 


*  i  i  i 

j=1  k=1  i=1  i  j 


In  (7.11)  Q 


(k) 


is  defined  by  (3.9). 


Now  we  claim  that  the  conditions  (7.8)  are  satisfied  for  f  given  by 


(7.11),  with  v  satisfying  (2.4)  and  a 
the  nonlinear  terms 


(3) 


satisfying  (3.8),  if  and  only  if 


,,  V  jV1  V  (r  *<:>)*  <i>D0OD(j)  (k)  (j)  .  (k)  (i)D(j)D(k)  (j)  (k). 

(7.12)  l  l  I  (L  $  A  R  R  a  ag  +  L  ^  A  R  R  a  J  afl  } 

3=1  k=1  i=1  i  3  i  3 


have  zero  mean  over  S(0).  The  reason  for  this  is  as  follows.  The  first 


three  terms  in  (7.11)  proportional  to  v  sum  to  zero  by  (2.4).  The  terms 


/  )r  ) 

linear  in  a'  and  the  nonlinear  term  proportional  to 

(k)  (k)  1  3  r  ( k )  2  .  c  f  rt \  i  a 

a  a^  =  —  -g-g—  [a  ]  have  zero  average  over  S(a)  unless  9^  is 

k  k 


/  v  \ 

constant  on  S(o).  If  0^  is  constant  on  S(0)  then  we  must  have  L  =  1/  , 
and  then  the  scalar  produce  of  the  sum  of  these  terms  with  L  is  zero  by 
(3.8). 

Thus  a  necessary  and  sufficient  condition  for  the  solution  derived  in 
section  3  to  satisfy  the  solvability  conditions  (7.8)  is  that  the  terms  (7.12) 
have  zero  mean  over  S(o).  However,  suppose  that  9^  and  0^  are  not 
constant  or  functionally  dependent  on  each  other  on  S(0).  Then  we  can  use 
=  an<*  ^2  =  as  tW°  coor<iinates  on  S(0),  and  the  mean  of  (7.12) 
over  S( 0)  is  clearly  zero. 

If  8^  and  8^  are  constant  or  functionally  dependent  on  each  other 
then  we  must  have  for  some  c.(x),  c.{x)  and  (y  (x) , . . . ,y  (x) ) 


3 


1 


°/)  *  •  J,  V*  • 


(7.13) 

Using  (7.10) 

(7.14) 

Since  (7.14)  holds  for  all  0,  in  particular  it  follows  that 


m  v 

c,8.  +  c6  =  y  y  y .  x  0 

3  3  fc  k  ^  1  l  lr  r 


(7.15) 


ci  =  l  , 

3  i=  i  * 


“tV  • 


HI  .  >  .  .  . 

But  by  the  definition  of  X  ,  \  X  4>  AU  has  left  null  vector  L  for 

rr  r=1  *r  x^ 

l  =  1,...,v.  Therefore  if  p  is  defined  by  (2.15)  using  (7.15)  the  matrix 

?  ( i) 

I  p  A  has  left  null  vector  L,  and  (2.16)  must  hold. 
i=1 
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This  proves  that  if  (2.15)  and  (2.16)  do  not  hold  at  any  x  for 

distinct  j  and  k,  then  the  mean  of  all  terms  (7.12)  over  S(o)  is  zero, 

and  v  and  a^^  satisfy  the  solvability  conditions  (7.8).  On  the  other 

hand,  suppose  (2.15)  and  (2.16)  do  hold  for  some  pairs  {j,k}.  Then  the 

solution  described  in  section  2  is  valid  if  and  only  if  for  each  pair  {j,k} 

the  average  of  the  term  (7.12)  is  zero  over  the  space  S (0)  corresponding  to 

the  left  null  vector  of  L  of  p.A  . 

i=1  1 

In  section  3  we  took  a  particular  form  of  solution  for  and  v*0^  to 
(3.3).  We  have  shown  that  this  is  sufficient  provided  that  re  ^  e  does  not 
occur.  To  treat  problems  involving  resonant  interactions  a  more  general  form 
of  solution  to  (3.3)  for  $  and  v^  must  be  used. 
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8.  Multiple  Characteristics 

In  this  section  we  generalize  the  results  of  section  2  to  the  case  of 

multiple  characteristics  for  a  large  class  of  systems.  This  class  includes 

systems  derived  from  a  set  of  conservation  laws. 

Let  us  regard  the  eiconal  equation  (2.6)  as  an  equation  for  $  ,  where 

X1 

we  can  suppose  that  is  a  timelike  coordinate.  Equation  (2.6)  is  an  m-th 

degree  polynomial  in  $  .  We  denote  its  m  real  roots  by 

X1 

g  (x,$  ,...,$  ).  Suppose  that  (2.1)  has  multiple  characteristics.  Then  v 

j  x2  xn 

of  the  gj's  will  be  the  same.  We  suppose  g^  =  g^  for  j  ■  1,...,v. 
Therefore  we  shall  take  $^(x)  *  $^(x)  and  0_.  =  0,^  for  j  =  1,...,v. 

We  assume  that  the  multiplicity  v  does  not  change,  and  we  also  assume  that 
there  is  no  resonance  i.e.  that  (2.15)  never  holds. 

We  use  (3.2)  in  (3.1)  and  use  the  result  in  (2.1).  Equating  to  zero  the 
coefficient  of  e°  we  obtain  (3.3).  Setting  each  term  in  the  sum  over  m 
equal  to  zero  gives  (3.5).  We  take  as  solutions  to  the  eiconal  equation 
(2.6),  ♦^1^(x),  $  ^V+1  ^  (x) , . . .  (x) .  To  simplify  notation  we  suppose  that 

the  other  families  of  characteristics  are  simple.  This  does  not  affect  our 
final  results  which  are  the  same  whether  or  not  there  are  other  multiple 
characteristics.  In  fact,  several  of  the  $^'s  can  correspond  to  the  same 
multiple  characteristic. 

Once  the  <J»^  are  determined  v^  satisfies  (3.5).  The  matrix 

l  a1  has  a  null  space  of  dimension  v.  We  denote  by 

i-1  xi  u 

{R^ 1 ^ , . . . ,R^V a  basis  of  the  space  of  right  null  vectors,  and  by 
{L(1) . L<V)}  a  basis  of  the  space  of  left  null  vectors.  We  write  v^0^ 


(8.1)  v(0)(x,8)  =  v(x)  +  l  a( j)(x,01)R< j)  +  £  a( j) (x,0 . )R( j)  . 


j*v+1 


Then  vv  '  clearly  satisfies  (3.5). 


The  only  difference  between  (8.1)  and  (3*6)  is  that  in  (8.1)  the  first 
v  a(j)  all  depend  on  the  same  fast  variable  8  . 

Equating  the  coefficient  of  e  to  zero  we  obtain 


l  l 


I  a11’?  +a‘1W(>'> 

U  V  II  V  11 


k=1  i=1  i 


(8.2) 


+  I  l  {(a(k)A(i)R(k)  ♦  a<k)^k)A(i)vR(k))  ♦  a(k>Q(k)> 


k-1  k-1 


y  n 

9k  -i  u 


*  l  I  l 

k-1  A-1  i  i-1  i 


In  (8.2)  8  is  taken  equal  to  8^  for  j  *  1,...,v.  now  equation  (2.4)  for 

v(x)  follows  from  (8.2)  by  averaging  it  over  9  provided  that  the  mean  of 
the  nonlinear  terms  proportional  to  a^aik^  is  zero. 

j 

In  fact  this  is  not  always  the  case  when  there  are  multiple 

characteristics.  The  mean  of  a^  (x,8^  )»gk^  (x^ )  is  not  necessarily  zero 

when  j  /  k,  and  1  <  j,  k  <  v.  in  general  the  terms  [a^agk*]av0  couple 

_  ^  ( i ) 

together  the  equations  for  the  mean  v(x)  and  the  amplitudes  a  J  (x, 61 ) . 

However  we  shall  show  that  the  mean  of  the  nonlinear  terms  is  zero  for 
systems  which  satisfy  the  following  condition. 

Condition  8.1.  For  all  1  4  j,  k  4  V: 


(8.3) 


I  ♦(1)A(i)R(j)R(k)  =  l  ♦(1)A(i)R(k)R(^) 
i-1  Xi  U  i-1  Xi  U 


If  there  are  other  multiple  characteristics  they  also  satisfy  the  analagous 
equation  to  (8.3). 

From  (4.14)  condition  8.1  is  always  satisfied  by  any  system  which  is 


derived  from  a  set  of  conservation  laws.  We  now  show  that  the  mean  of  the 
terms  a^^agk^  is  zero.  Trouble  can  only  arise  when  1  4  j,  k  4  v  so  that 
both  a  ^  and  a  ^  depend  on  0  . 


If  j  =  k  then 


Thus  a^^ai^  is  an  exact  derivative  and  has  zero  mean,  since  a  ^ 


is  a 


bounded  function  of  0  .  If  j  /  k  we  consider  the  terms  proportional  to 


a<*>a<k)  and  a(k)a^> 
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From  (8.2)  these  are 


,Q  ..  v  . (1)  (i)  ( j)  (k)  (j)  (k)  A  r  .  ( 1 )  (i)  (k)  ( j )  < k >  ( j ) 

(8.4)  I  <|i  A  R  J  R  a  aA  +  l  *  A  R  R  J  a  aflJ 


i-1  Xi  U 


.  ,  x.  u 

i=1  l 


Using  (8.3)  we  rewrite  (8.4)  as  an  exact  derivative 


(8.5)  i  3§-  (.(3>.ik>) 


.  4  x.  u 
1=1  1 


Therefore  (8.4)  averages  to  zero.  This  shows  that  v(x)  satisfies  (2.4). 


To  obtain  equations  satisfied  by  a 


(j) 


for  j  *  1 , 


,v  we  average 


(8.2)  over  ( ®v+i • • • • #®m)  and  take  the  scalar  product  of  the  resulting 
equation  with  L^P^  (p  *  1,...,v).  Then  we  find  the  same  equations  as  those 
obtained  by  Choquet-Bruhat  [  1  ]  : 


(8.6) 


l  l  L(p)A(i)R(j)a(j)  +  l  l  l  L(p)*M)A(i)R(k)R(j)a(k)a<j) 
j=1  i=1  Xi  j=1  k=1  i=1  xi  U 


v  n 


+  l  l  L(pV1,A(i)vR(3)a‘3)  +  l  L(p)Q(j)a(j)  =  0  . 


j=1  i=1 


x,  u 

l 


j=1 


We  have  the  following  lemma: 

Lemma  8.1.  There  are  scalar  functions  c^x),  for  i  =  1,...,n,  such  that 
for  all  j  =  1,...,v  and  p  =  1,...,V; 

L(p)A(i)R(j)  =  ciL(p)R<;j)  . 

The  proof  of  lemma  8.1  follows  that  of  Lewis  [12]  almost  exactly. 


Using  lemma  8.1  in  (8.6)  we  find 


(8.7)  l  l  L(p)R(j,a^)  *  l  l  l  l(*V1V1,R(k,R<J,a<k,.W 
j=1  i-1  0  j=1  k-1  i-1  Xi  U  91 


v  n 


+  I  I  ♦  I  l<p>q<3).(J>  =  o  . 


j=1  i-1 


X.  U 


j-1 


(p  —  1,...,v)  . 


In  equation  (8.7),  0  is  defined  by 


(8.8) 


^i 

IT  “  ci(x)  * 


Equation  (8.7)  is  a  system  of  quasi-linear  equations  for  the  v 
scalars  a^)  in  two  independent  variables  a  and  0^. 

Equation  (8.7)  reduces  to  a  system  of  ordinary  differential  equations  in 
the  special  case  that  there  are  A(x)  e  R  and  y(x)  e  R  such  that 


I  iW'V^rWrW  -  AlW”*  , 

i-1  xi  u  Jx 


(8.9) 


l  . 


i-1  xi  u 
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Appendix 

In  this  appendix  we  specialize  the  results  of  section  2  to  the  gas 
dynamics  equations 


(A.  1 ) 


+  div(Pu)  =  0 


pufc  +  Pu*Vu  +  Vp  =  0 


p  =  <P 

Equations  (A.1)  are  the  equations  of  motion  for  the  isentropic  flow  of  a 
compressible,  inviscid,  ideal  gas.  In  (A.1)  p  is  the  gas  density,  p  the 
pressure  and  u  the  velocity  of  the  gas.  The  quantities  x  and  are 
constants . 

The  weakly  nonlinear  geometrical  acoustics  solution  (2.2)  to  (A.1)  for 
x  and  u  in  R3  is 


(A. 2) 


■> 

p 

P  (x,t) 

P(x,t) 

= 

+  e 

u 

U(x,t) 

u(x,t) 

m 


+  e  l  a(3)(x,t,<|>(3)  (e)R(3>(x,t) 
j=1 


m* 


+  £  i  {a(k'1)(x,t,i|;(k)/e)R<k'1,(x,t)  +  a(k,2)  (x,t,*(k)/e  )R(k'2 }  (x,t ) } 
k=1  2 

+  0(e)  . 

In  (A. 2)  (Pg ,U)  is  any  exact  solution  to  (A.1).  The  mean  value  of  the 
perturbation  in  (P,u)  about  (pQ,U)  is  (ep,eu).  The  sum  over  j 
represents  a  superposition  of  high  frequency  sound  waves,  and  the  sum  over 
k  is  a  superposition  of  high  frequency  vorticity  waves. 

The  means  (P,u)  satisfy  (2.4)  which  gives  the  acoustics  equations 


(A. 3) 


pfc  +  div(PQu  +  PU)  =  0 


Vt  +  v«#>  +  P°T  +  P0«*?u 
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In  (A. 3)  Cq  is  the  linearized  sound  speed 
(A. 4) 


=0-nr’)V2  . 


and  is  the  derivative  taken  along  the  streamlines  of  U 

(A. 5)  |_=|_+U.V  . 

The  eiconal  equation  (2.6)  for  (A.1)  when  x  e  R°  is 
(A. 6)  (♦£  -  cq|V<M2)^_1  =  0  . 

The  phase  functions  and  in  (A. 2)  satisfy 


(A. 7) 


(*Tj>)2  "  co|7*<j)|2  = 


We  consider  the  sound  waves  (which  depend  on  $  J  )  and  the  vorticity  waves 


- » 
in 


(k) 

(which  depend  on  <|i  )  separately. 


For  the  sound  waves  the  rays  (2.7)  corresponding  to 


,<J> 


are  given  by 


(A. 8) 


dt 
ds . 


-  u 

j  T  dS j 


-  Vj43>u  -  »Jjc^(j)  , 


and  the  vector  R^  is  given  in  terms  of  $ 


(j) 


by 


(A. 9) 


R^>  =  X  . 

3 


-p  c~Vj) 

0  0  x 


7* 


(j) 


The  scalars  Uj(x,t)  and  X^(x,t)  in  (A. 8)  and  (A. 9)  are  arbitrary  functions 
which  may  be  chosen  in  whatever  way  is  convenient. 

The  amplitude  a^  ^  (x,t,$^  Ve)  of  the  sound  wave  is  given  by 

(2.9).  For  (A.1)  Ejfs^Bj)  is 


(A. 10)  E 


(j) 


exp 


-L3  i 


TT 


2  -  ^  43>diV  0  "  2^  VP0*7*<3>},JjdSj 


The  integral  in  (A. 10)  is  taken  along  the  rays  (A.8).  The  first  term  in  the 
integrand  gives  the  change  in  the  wave  amplitude  due  to  changes  in  the  ray 
geometry.  The  second  two  terms  in  the  integrand  give  the  change  in  the  wave 
amplitude  due  to  inhomogeneities  in  the  zeroth  order  flow  p  =  p^,  u  =  U. 
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The  modified  phase  function  in  (2.9)  is  defined  implicitly  by  (2.11) 


which  becomes 


s  . 


=  d>(j)  -  e  1=1  F.(@.,C./e>  /  j  X.^Vc'Vli.ds’. 

2  3  3  3  J  sQ  j  T  0333 


(A. 11) 


-  e  /s^  ^j)u‘V^(j)lJjds'  +  e  —  /s^i:'))2"P01^dsj 


If  shocks  form  they  are  fitted  into  the  solution  using  (2.13). 

(k) 

The  phase  function  ip  for  the  vorticity  waves  corresponds  to  a 
characteristic  of  multiplicity  n  -  1  when  x  e  R°.  For  three-dimensional 
flows  the  characteristic  has  multiplicity  two.  Therefore  as  explained  in 
section  8  we  obtain  two  null  vectors  R^k/1'  and  R(k,2)  for  and  two 

amplitudes  a^k<1^  and  a^k'2^  which  are  coupled  together.  The  equations 
satisfied  by  a(k,1)  and  a(k'2)  happen  to  be  linear  for  (A.1). 


The  null  vectors  R^kf1^  and  R^k'2*  are  given  in  terms  of  ip 


(k) 


by 


(A. 12) 


R(k'1)  -X. 


0 

0 

'I' 

*2 

0 

-*<k) 

*<k> 

X1 

X3 

0 

Vk) 

-  — 

-  X2. 

In  (A. 12)  X^(x,t)  is  an  arbitrary  scalar.  To  put  the  equations  for 

( k  1 )  ( k  2  ) 

a  '  and  a'  '  in  their  simplest  form  it  is  convenient  to  let 


(A. 13) 


(k) 

W 

1 

\<k> 
T  v 

0 

2 

(k) 

W2 

=  X  a(k'1) 
Aka 

_Jk> 

X1 

♦  Xa(k'2) 
k 

X3 

(k) 

w 

3 

0 

w- 

-*<k) 
L  X2 
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and 


Then  w^k'  (x,t,^k  Ve )  is  a  velocity  vector  orthogonal  to 
satisfies 

f  V  J  (  v ) 

(A., 4)  »‘kl  *  4  v<k).Vu  -  *!"  •7"  V*(k)  -  0  . 

T  0  |V*lKY 

(k )  ( k ) 

We  solve  (A. 14)  for  w  (x,t,o)  and  evaluate  a  at  ^  /e.  when  U  is 

(k)  (k) 

independent  of  x  (A. 14)  gives  three  uncoupled  equations  for  w.^  ,  v#2  and 

(k) 


We  may  want  to  use  the  local  sound  speed  c(x,t)  as  a  dependent  variable 
instead  of  p(x,t).  They  are  related  by 
(A. 15) 

Then  we  seek  a  solution 


2  Y-1 

c  =  <YP 


(A. 16) 


f 

c 

'cQ(x,t) 

c(x,t) 

= 

+  e 

_U(x,t)  j 

.u(x,t) 

m 


+  e  l  a(j>(x,t,<J>(:i)/e)R(  j)(x,t) 

j-l 


m 


+  E  I  {a(k'1)(x,tyk)A)R(k'1)  +  a(k'2)(x,t,«|)(k,/e)R(k'2)}  +  0(e2,  . 

k=1 


In  (A. 16)  c  =  Cq,  u  =  U  is  an  exact  solution  of  the  gas  dynamics 
equations.  The  means  c  and  u  satisfy 


(A. 17) 


_  v«i  —  v_ ^  _  . 

c  +  -  c  div  u  +  -  c  div  U  +  u*^c_  =0  , 

T  2  0  2  0 


uT  +  j-  V(c0c)  +  u*Vu  =  0  . 


The  R^^  in  (A.  16)  are  given  by 


(A. 18) 


R(j)  =  X. 


-(T-1)43>/2c0 


(j) 


where  A_.(x,t)  is  an  arbitrary  scalar.  The  remaining  equations  are  exactly 
as  before:  satisfies  (A. 7),  a^  is  found  from  (2.9),  (A. 10)  and 
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(A. 11)  with  P/Pq  ■  c/Cg,  satisfies  (A. 7),  rO^#!)  an<j  r^'2^  are 

given  by  (A.  12)  and  w^k*  satisfies  (A.  14). 

The  equations  for  plane  flow  (x  e  H2)  are  obtained  from  those  above  by 
suppressing  the  x-j-dependence  and  setting  a*  '  '  «  0.  The  equations  for  one 
dimensional  flow  follow  by  suppressing  the  X2  and  x^-dependence  and 
setting  a(k,1)  =  a(k,2)  «  0. 
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